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Abstract 

We formulate a general theory of the dielectric response of a lattice with a 
structural transition and a polarisation instability due to a soft-optic mode 
coupled to low energy electronic excitations by the electromagnetic fields. 
The electronic excitations considered are in two-limits; those of a low density 
of free-electrons or those of a low density of strongly localised electrons in 
the Coulomb-glass phase. The ferroelectric transition in the absence of the 
electronic-excitations and the low energy dielectric properties are shown to 
be strongly modified. 
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I. INTRODUCTION 



In the last few years there has been an enormous revival in the study of ferroelectric and 
ferroelectric materials 0, 0. In a large part it is due to the investigation of complex 
materials with very large number of atoms per unit cell and substantial disorder. The 
disorder also introduces low-energy electronic states. Such materials display what has come 
to be known as 'relaxor-ferroelectric' properties. This term appears to cover all properties, 
not seen in simple ferroelectrics- for example BaTiO^. Generally these materials do not 
have a ferroelectric transition but do have a very large dielectric-constant, the polarization 
displays complicated hysteresis loops and has a strong frequency and history dependence at 
low energies. 

Various ideas for the form of disorder and its effect on the ferroelectric transition have 
been proposed 0], 0, 0. But to our knowledge, the important effect of low-energy particle- 
hole excitations due to shallow electronic levels has not been treated. Ferroelectricity (like 
superconductivity) is rather special because the order parameter, the polarization, couples 
to the low-energy electronic excitation through the electromagnetic-fields. In this paper 
we study such effects. The effects of disorder are neglected here by adopting a virtual- 
crystal aproximation. In fact both effects, disorder and low-energy electronic excitations are 
important; we hope to treat them together in the future. 

This paper is organised as follows. In Section 2, we present a general Lagrangian for a 
polarization instability coupled to electronic excitations through the electromagnetic fields. 
The generalization to the disordered case is also indicated. In section 3, we consider the 
simple case that the low energy excitations are due to a low density of free electrons. While 
it is common-place knowledge that free-electrons will screen the electric field so that macro- 
scopic fields cannot arise, it is worthwhile to precisely derive the conditions for this and 
see the condition for the formation of ferroelectric domains at low enough charge density. 
The domain walls in this case may be charged. We also derive the low energy charge and 
current response functions for this case. In Sec. 4 we take up the more interesting problem 
of localised electrons in a Coulomb glass phase and derive the charge response functions. 



II. THE LAGRANGIAN FOR SOFT OPTIC MODES IN PRESENCE OF 
LOW-ENERGY ELECTRONIC EXCITATIONS 



We consider the following Lagrangian density (at imaginary times), L = Lferro + Lem + 
Lint + where 



-'ferro 
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is a Langrangian describing the polarization vector P. For uiq = a{T — Tc) Eq. (la) describes 
an instability of the system (in the absence of coupling to electric fields) when T — > Tc. 
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is the usual Lagrangian for EM fields in terms of vector potential A and scalar potential 
(f) in the Coulomb gauge, V.A = 0. In the microscopic theory, ujq is the frequency of the 
zone-center transverse optic phonon. The polarization vector P couples to the electric field 
E through the usual Lagrangian Lint, where 



-'int 



P.E = -P. 
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In the presence of charges, the EM potentials pick up additional terms in their La- 
grangian. In the Coulomb gauge they are 



= -- 
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where xl and xt are the longitudinal and transverse response functions of the electrons to 
scalar and vector potentials and A, respectively. In Eqs. (1), P, and A depend on the 
position r. For the translational invariant case, this may be directly Fourier-transformed. 
In the case of disorder, they are random functions of position and the response functions 
Xl and xt are functions of two co-ordinates. In what follows, we shall assume an effective 
translational invariant approximation, where disorder is considered averaged over at the 
outset. 

In the Coulomb gauge, the Lagrangian can be separated into longitudinal and transverse 
parts, L = Ll + Lt + Lu- We obtain after fourier transforming. 



(2a) 



where v{q) = Aix/q is the Coulomb potential, and 
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AojA, 
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where we have divided the polarization vector P into longitudinal and transverse parts 
P = Pl + Pt, and 



u 



L^ = -{Pl + Pt) 
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(2c) 



is the fourth-order term in Lferro- Notice that L„ mixes the longitudinal and transverse 
parts of P. 

To study the ferroelectric instability, we take the iu = Q limit and first minimize the free 
energy with respect to the and A fields. We obtain 



iq.Pi 



iuoPx 



(3) 



A ^ . 

c {iuY — c^q^ + AnxT^q,^) 



0. 



Substituting these back into Eq. (|^), we obtain an effective free energy in terms of P only, 
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In an infinite medium Pt decouples from EM field completely in the static limit, as is 

— * — * ^ 

required by Faraday's Law that B — V x E when a transverse electric field exists. This 
is not true in a finite sample, where a macroscopic Pt, if it exists is discontinuous at the 

— * 

surface. A surface charge - density a — P.n is then generated leading to an electric field. 



III. FREE ELECTRONS 

For free electrons, we may use the fact that 

ivi^r - XLmr ^ (5) 

in the small q limit. In Eq.(5), Qsc is the static screening length, given in terms of density 
by the Thomas-Fermi length at low temperatures and by the Debye length in the classical 
regime at high temperatures, i.e. 

qI ~ Qtf = QTrnoeyEf = 4(3/7r)V3nyV«o (6) 

for T << Ef. Here Uq is the density of electrons, Ef their Fermi-energy and Oq the Bohr 
radius, and 

Qsc ~ ?D = 47rnoeVA;r. (7) 

at high temperature. For most problems of interest in ferroelectrics, the density is low 
enough that near room temperature the Debye Screening is applicable. 

In the absence of free electrons, an instability developes in Pt field in the q ^ limit as 
T Tc (or ujo — ^ 0) always before any possible instability in Pl. This instability occurs at 
Tc, where the transverse optic phonon frequency — > 0. The magnitude of Tc will in general 
change in the presence of free electrons due to microscopic effects not treated here. The 

— * — * 

presence of free electrons bring the transition in Pt and P^ to the same temperature Tc 
because of screening, i.e. in infinite samples, no uniform electric-field can be produced due 
to a macroscopic P because the electric field sets up by the ions oscillation is screened. The 
electric field in the ordered state of the system with < is given by 

E = -V < 6 >= Urn' ' 



<z^o yv{q) 1 - XL J ' 

Here a is the surface charge density at the surface of the sample at which Pt has a normal 

— * — * — * 

component Pt-u ^ 0. E goes to zero as long as 0) is finite, i.e. < Pt >^ does 
not lead to uniform electric field in the presence of free electrons. This does not necessarily 
mean that ferroelectricity is destroyed by arbitrarily small density of electrons. As shown 
below a domain structure of electric-fields is still possible provided the charge density is low 
enough. 
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At temperatures T < Tc, it is easy to minimize the energy to show that < P >— Pg — 



2 — * 

The Langrangian can be written as L = L{Po) + L{P'), where P' are fluctuations 



around P = P^. We find that to second order in P', L(P') = Ll(P[) + Lt{Pt), where 



Ll{P'l) = 



X 



r + Ts 
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and 



LriP^) = —. 
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(8b) 



Notice that the only difference between T > Tc and T < Tc in the quadratic fluctuations in 
P is that uJl is replaced by — 2a;^ as temperature decreases from T > Tc to T <Tc. 

IV. DOMAIN WALL FORMATION 



For usual (insulating) ferroelectrics with finite dimensions, domains are formed to min- 
imize the electric field energy ~ / d^x^, coming from effective surface charges a[x) ~ 
n.P{x), where n is a surface unit vector. In this section the criteria of domain formation in 
the presence of free charges is studied. For simplicity we shall consider a ferroelectric with 
a finite width 2L in 2;-direction, and with infinite extent in x — y plane. We shall assume a 

— * — * 

transverse polarization field Pt with Pt pointing in ^-direction, and with magnitude func- 
tion of X— and y— coordinates, i.e. P = zP{x,y) for —L < z < L. For a finite L surface 
charges o-{x,y) = {6{z — L) — 6{z + L))P{x,y) will be set up, and the total energy of the 
system is given in momentum space by 



(9) 
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+^f^dz I I dxdyP'', 



where a{q) = {2i)sin{qzL)P{qx.,qy) is the fourier transform of the surface charge density. 
Notice that q.Priq) = in the sample. Minimizing the energy with respect to 4> we obtain 



l {2sin{q,L) \P{q^.qy)\f ^ ,2 , ^2^2n| p. 



(10) 
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-\ f_^dx j j dxdyP\ 

Summing over g^, and using the results that 
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ql + Q^ \/-oo^^^2(g2 + g2) 2g' 
when L — > oo, where is a positive number, we obtain 
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(11) 
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dxdy-^P^, 



where q = (qx^Qy) in Eq. (|TT]). Minimizing the free energy with respect to P{q) and q we 
see that the system has now an instabihty at finite q which indicates the instabihty of the 
system to domain formation. Expanding the quadratic terms in P in Eq. (|ll]) in powers of 
at small q, we find that the system has an instability towards domain formation when 



11 
X 



2tt 

qlcL 



< 0, 



(12) 



or when ^ Lj^ = (2^t^)^- It is easy to show from Eq. (pH]) that in the absence of 
free electrons (g^c = 0) gives the characterestic size of the domains. The criteria for 
formation of domains (0) can thus be interpreted as simply saying that the screening length 
of free electrons should be longer than the natural domain size in order for domains to be 
formed when free electrons are present. 

The presence of a finite charge density has important consequences for the switching 
and hysteresis properties of domains by reversing external-fields. Domain reversal must be 
accompanied now by motion of charges. The time scale for this is now limited by the time 
for the charges to diffuse over the characterestic size of the domains. 



V. DIELECTRIC FUNCTIONS 

To obtain the dielectric function at T > Tc, we neglect the Lu term and integrate out the 
P fields first in the Lagrangian (^. As a result, we obtain two effective Lagrangain Lf.ff{(j)) 
and Leff{A) for the electromagnetic potentials. The dielectric functions can be obtained by 
identifying 

LefM = -^^4> , (13) 

r fX\ eT(g,i^)[(^^)'-cV] ^2 

^eff[A} - — A^, 

where we have again considered the Coulomb gauge V.A = 0. After some algebra we obtain 

4:7rx 

eL{q,iuj) = l-v{q)xL{q,i^) - j-^ TTT^ (14a) 



and 



eriq, lu) = 1 — — 14b 
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where uJo{(l)'^ = + 7^0'^- ^ similar expression can also be derived for T < Tc. The only 
difference is that uj^ is replaced by —2u!l when T < T^. Notice that uJo{qY is positive definite 
at all temperatures T 7^ T^. 

The longitudinal electromagnetic collective mode dispersions in the system can be ob- 
tained by solving the equation 6^(5, co'(g)) = 0. In the small q limit and with uj » q, 

v{q)xL{<l,^^) TjT, where up is the electronic plasma frequency, and we obtain the secular 
equation, 

^ 2 2 Tyl = ^- \^^> 

(This equation is not modified in the classical limit for long- wavelengths.) Solving the 
equation we find that there are in general two solutions ui and uj2 in the equation. At small 
q and for uJoiqY « ujp, we obtain 

iof -^ujI + Attx (16) 
\ujf, + A7TxJ 

The transverse electromagnetic collective mode dispersions are obtained by solving 
uj'^eT{q,uj) — (?q^ = 0. We obtain after some simple algebra a secular equation very similar 
to Eq. ([15|), except that Up is replaced by Up + c^q"^, with similar solutions 



^1 



2 .^| + 47rx + cV (17) 

^2 ^ 2 I 2 2 1/1 • 

Notice that uJi gives the usual plasmon and polariton modes, respectively for longitudinal and 
transverse EM fields in metals, with modifications coming from polarization field, whereas 
UJ2 gives the collective modes in longitudinal and transverse polarization fields Pl and Pt, 
respectively, with modifications coming from metallic behaviour. Notice that at small q and 
when — > 0, ci;2 becomes very small and will be in the hydrodynamic regime where xl 
would have a diffusive form 

dn Dq^ 



dfx iuo — Dq^ 

In this regime, the dispersion uj2 will be modified. In particular, a finite damping of the 
collective mode will occur because of coupling to the electronic degrees of freedom. 



VI. FERROELECTRICS WITH LOCALIZED ELECTRONS 

We shall consider here ferroelectrics in the presence of strongly localized electrons in 
the Coulomb glass phase. In the dilute particle-hole pairs approximation (see Appendix), 
the localized electrons have only longitudinal response to external electric fields, and with 
response functions XcciQi^) given in the Appendix. Their effects on the ferroelectrics can 
be studied by replacing xl for free electrons by Xcg- First we consider = and study 
ferroelectric instability and domain formation at different temperature ranges. 
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A. ferroelectric instability and domain formation 



As a function of temperature, the response function XcdQi 0) has three different regimes. 
First we consider high temperature ksT >> A where the Coulomb glass effect is unimpor- 
tant and the response of the electrons are dominated by phonon-assisted hopping. This 
is probably the regime of experimental interests in ferroelectrics. Using Eq. ( |A11|) for 



R^Xcg{(1i^) in Eq.(lO) and repeating the analysis afterward we find that the free en- 



ergy is given by Eq.(ll) with \Jq^ + q1c replaced by q^/e in first term of Eq.(ll), where 

e ~ eo + C^"^i44o"'° )(fc^)^) where eg includes all the non-singular contributions to e. The 
characteristic domain size L£, is given by 

~ (?2!^)i, (18) 

which differs from the pure ferroelectric case only by the factor e. Notice that e is expected 
to be large in the Coulomb or fermi glass regime and we expect that the effect of introducing 
localized electrons in the system is to enlarge the domain size. Notice that for small enough 
L, we may enter a regime where L^jao « (rT/ao)"^. In this regime, we should use Eq. 
( |A1CI| ) for Rexccil) 0). In this regime, we find that the enhancement in domain size by the 
factor e is lost, and the analysis is similar to the case of metals except that the screening 
length qj^ is given by 



144 



Notice that in contrast to metals, the screening length increases as temperature is lowered. 

As temperature is lowered we find that for large enough L the domain size is still given 
by Eq. (^) except that the value of e goes thought two more regimes. At temperature 
A > ksT but ujt > I(ro), we find 



whereas at ujt < H^o), we find 



2 4 4 ^-1 

e ~ eo + 1^) ' 



To 

e ~ Co + ——ln{^). 



IOtt ' T 



Notice that our theoretical prediction of domain size variation as a function of temperature 
can be tested experimentally. Note that for small enough L, the effect of screening will 
become stronger as in the high temperature case fc^T > A, and the enhancement factor e 
in domain wall size is lost. 

It should also be kept in mind that the Coulomb glass has hysteretic low frequency 
response. One should therefore expect a slow and sluggish transition over a wide temperature 
range starting from about the to roughly when uJq + Attx ~ 0. 
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B. (longitudinal) collective modes 



At T > Tc, the longitudinal collective modes of the system can be obtained by solving 

1 - v{q)xcG - 2 ^""^i ,2 = 0- (19) 
In the >> g limit, we obtain the secular equation 

= ^o{qf + (20) 

where Xccil 0,0;) are given by Eqs.(A6a), (A6b), at temperature range /cbT < A at 
u < I{ro) and u > I{ro), respectively. At temperatures ksT » A, Xccil ~^ 0,0) is given 
by Eq.(A9). Although the quantitative behaviours of the collective modes are different, 
logarithmic corrections appear in the pure dielectric collective mode frequency uJo{q) in all 
the above cases. The logarithmic correction vanishes in the ferroelectric instability limit 
u ~ ci;o(0) — >■ 0, in agreement with the static instability analysis which indicates that the 
ferroelectric instability occurs exactly when ci;o(0) — 0. Notice, however that a collective 
mode solution with cu'^ ~ ^oiq^ + ^ exists even up to the instability limit cOoiO) 0. The 
polarization collective modes are in general much broadened by the Coulomb glass particle- 
hole excitation spectrum because of the very slow vanishing of ImxcG{<li^) at small uj. 
Notice that at T < we may obtain the collective modes by replacing ujI by —2uj1 in Eq. 
(p!9|). This is valid as long as we are interested at responses of the system at wave vectors 
q < L^K 



APPENDIX A: ELECTROMAGNETIC RESPONSE OF COULOMB GLASS 

In this appendix we summarize our theoretical understanding for response functions of 
Coulomb glass systems. (Notice that we have set = 1 in the main text. Here is 
kept to make the equations more clear). The electromagnetic responses of Coulomb glass 
systems are rather complicated, and can be divided into many different regimes, depending 
on temperature, frequency and wavevector. A detailed discussion of these can be found in 
ref[ i]. 

The starting model for responses of Coulomb glass systems is the effective two-site Hamil- 
tonian by Shklovskii and Efros [0]. The Hamiltonian is derived under the assumption that 
the system is in the strongly disordered regime where the electrons are strongly localized. In 
this regime tunneling events of electrons from one site to another is rare and can be treated 
in a dilute particle-hole pairs approximation where interaction between different tunneling 
events are neglected and the tunneling of electron between two sites (particle-hole pair) is 
described by a two-site Hamiltonian 

Hij = EiUi + EjUj H Tiirij + I{rij){afaj + a'j'ai), (Ala) 

where rii is the electron occupation number (= 0, 1) at site i, rij is the distance between the 
two sites, /(r) ~ Ioe~'^/°'° is the effective tunneling matrix element between the two sites. 
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where Qq ~ lattice spacing and 0(0"*") are electron annihilation and creation operators, the 
on-site energy Ej, is given by 

E^ = e,+ Y: T^^, (Alb) 



where —A < Si < A is the random on-site energy with bandwidth 2A and the second term 
comes from Coulomb interaction. Notice that electron spins are ignored in the theory. Hij 
can be diagonalized easily to obtain the following energy spectrum, 

Et^l{E, + Ej)±^ (A2) 

E2 — Ei + Ej H 

where E^ and E'f are the two possible energy states when the two-site system is oc- 
cupied by a single electron, and E2 is the energy when both sites are occupied, F = 



Ei — EjY + AI{rijY. Notice that for the pair of sites to be occupied by a single elec- 
tron and participate in the electromagnetic response, we must have E^ < ^ < E2 — E^ , 
where n is the chemical potential. The whole system is described as a diluted gas of indepen- 
dent (i.e. non-interacting), randomly distributed particle-hole pairs and the total response 
of the system to external electromagnetic fields is given by by sum of the response of in- 
dividual pairs. For example, the density-density response function of a pair in the two-site 
Hamiltonian approximation is found to be (at temperature T — 0) 



Xp(g,^w) = 2— — -(1 - cos{q.rij)) 



(A3a) 



F(n,-)' ' [ico - r{nj) iuj + F(n,)_ 

and the average density-density response function of the whole system is 

XcG(g, iuj)^<Y^ Xp{q, i^) > (A3b) 
p 

rA rA r T ( rV 2T(r) 

where F{Ei, E2,r) is the probability that a particle-hole pair with on-site energies Ei and 
E2 is found separated by distance r. For Ei, E2 larger than the Coulomb gap energy A ~ 
(e^/ao)^/^/A^/^, F where g — n/Aal is the density of states above the Coulomb gap, n is 
the average number of electrons per site. For Ei, E2 much smaller than A, F ~ p{Ei)p{E2), 
where 

P{E) = —,E\ 

Notice that a two-site system behaves as a electric dipolc and has only longitudinal response 
to external electric fields. Evaluating the integrals, we find for I{ro) « uj « A, 

Imxcciq,'^) ~ ao7r5(^eV^(l - (A4a) 
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where = aoln{^), Tq = e^/A. For lu « I{ro), we obtain 

ImxcG[q,^) ~ TTT^^K^ )• A4b 

The real part of the response function can be obtained using the Kramers-Kronig relation, 
we obtain for qro << 1, 

ReXcGiq,0) ^ ^^^;^q'^ln{qro), (A5a) 

and for qro >> 1, 

Rexcciq, 0) ~ -^e'r^ = -yr^A. (A5b) 



For /(ro) > uj and gr^ << 1, we obtain 



and for uj > I{ro), 



ReXcGiq,^^) ~ ^5^9 ln{ln{—)), (A6a) 



ReXcG{q,uj) ^l'^^^- (^6^) 



At finite temperatures T < A, eq. ( |A4b ) for Imxcci^^^) is cutoff at a low frequency 

1 

~ e~^^°/^)' by the thermal mechanism of Mott- variable range hopping ^ (in the pres- 
ence of Coulomb gap) where To ~ (? jao- Correspondingly, at temperature ujt < -^(^o) 
ReXcG{q,0) is replaced by [§ 

«n-ca(,. 0) ~ ^.iM^) ~ (A7) 

at a; < Ut, {q = 0) and at qvo < Vo/r^^, (a; = 0). At temperature ojt > I{^o), a corresponding 
analysis gives 



Rexccio, 0) ~ -^gV^, ~ y)'' (A8) 



for the same low frequency range and at small wavevector range qvo < {to/t^^^Y'. 

At still higher temperatures fc^T >> A where the Coulomb gap becomes unimportant 

and phonon-assisted hopping becomes dominant, we enter another regime where xcg become 

1 

different again At frequencies qr^ « 1 and uj > Cot, where Cot ~ ujpe~^'^^''^'^^^ , where 
is a characteristic phonon frequency, we obtain 



RexcG{q.uj) [^^^ q'rl (A9) 



and for u = 0, 
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Rexcciq,^) - -^aliksT). (AlO) 

For frequencies oj < Ut and for wavevectors qao « (ao/rx)^, where ~ (fc^)^O'o is the 
characteristic thermal hopping distance, RexcG saturates at a value ^ 

ReXCG{q,0) (All) 

Notice that the Mott variable range hopping introduce characteristic cutoff length scales 
r^y) and at different temperature regimes. As a result the dielectric constant e{q ^ 
, 0) = 1 — v{q)xcG{<l 0, 0) becomes finite but large at low temperatures although the 
temperature dependence of its value goes through several different regimes 0. 
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